Constraints on the origin of the ultra-high energy cosmic-rays using cosmic diffuse 

neutrino flux limits: An analytical approach 



Shigeru YoshidaEI and Aya IshiharaQ 

Department of Physics, Graduate School of Science, Chiba University, Chiba 263-8522, Japan 

(Dated: February 17, 2012) 

Astrophysical neutrinos are expected to be produced in the interactions of ultra-high energy 
cosmic-rays with surrounding photons. The fluxes of the astrophysical neutrinos are highly depen- 
dent on the characteristics of the cosmic-ray sources, such as their cosmological distributions. We 
study possible constraints on the properties of cosmic-ray sources in a model-independent way using 
experimentally obtained diffuse neutrino flux above 100 PeV. The semi-analytic formula is derived 
to estimate the cosmogenic neutrino fluxes as functions of source evolution parameter and source 
extension in redshift. The obtained formula converts the upper- limits on the neutrino fluxes into 
the constraints on the cosmic-ray sources. It is found that the recently obtained upper-limit on the 
cosmogenic neutrinos by IceCube constrains the scenarios with strongly evolving ultra-high energy 
cosmic-ray sources, and the future limits from an 1 km 3 scale detector are able to further constrain 
the ultra-high energy cosmic-rays sources with evolutions comparable to the cosmic star formation 
rate. 

PACS numbers: 98.70.Sa, 95.85.Ry 



I. INTRODUCTION 

The origin of the ultra-high energy cosmic-rays (UHE- 
CRs) has been a long-standing important question in 
astrophysics. While the observations by Auger [l|, Q 
and HiRes Q indicate that cosmic-rays with energies 
above ~10 18 5 eV are of extra-galactic origin, identifica- 
tion of astronomical objects responsible for the UHECR 
emission has not been achieved. Neutrinos, secondary 
produced by UHECR nucleons, are expected to provide 
the direct information on the UHECR origin, since a 
neutrino penetrates over cosmological distance without 
being deflected by cosmic magnetic field nor absorbed 
by the photon field. In particular, the intensity of the 
"cosmogenic" neutrinos [J] produced by the collisions 
of UHECR nucleons with the cosmic microwave back- 
ground (CMB) photon via photo-produced 7r meson de- 
cay as 7r ± — > l^ v ii s^Ve^fj,, known as the Greisen- 
Zatsepin-Kuzmin (GZK) mechanism [5|, indicates red- 
shift distributions of the parent UHECR sources [f| 0] in 
the Universe. The source distributions derived from the 
cosmogenic neutrino intensities can then be compared 
with distributions of known classes of the astronomical 
objects possibly responsible for the UHECR emissions. 
Therefore reliable extractions of the UHECR source dis- 
tribution function (s.d.f.) is one of key issues in cos- 
mogenic diffusive neutrino searches. Constraints on the 
sources of UHECRs derived from the measurements or 
upper-limits of the ultra-high energy neutrino flux are 
highly complimentary to the constraints from the diffuse 
photon flux [1, Q , because the former does not rely on 
uncertain estimation of extra-galactic background light 
(EBL). 



In this work, we develop a method to bound the 
UHECR source evolution and its redshift dependence in 
a comprehensive way without introducing specific astro- 
nomical models. We derive an analytical formula to cal- 
culate intensities of the neutrinos produced by the GZK 
mechanism in the range between 100 PeV and 10 EeV. 
Using the formula, we extract the relation among the 
neutrino intensity and the UHECR s.d.f. parameters. 
The use of the analytical formula allows us to calculate 
neutrino intensities in the full phase space of the source 
evolution parameters without an intensive computational 
task. The analytical formula can also be used as a practi- 
cal tool to approximately calculate cosmogenic GZK neu- 
trino intensity with given UHECR s.d.f., for example, for 
the performance studies of the future detectors such as 
KM3NET [10]. Finally we present model independent 
constraints on the UHECR sources using the obtained 
formula with the upper-limit [ll| and the future sensi- 
tivity [12] on the cosmogenic neutrino detection by the 
IceCube neutrino observatory [l3j . 

The standard cosmology with H ~ 73.5 km sec -1 
Mpc" 1 , Q M = 0.3, and Q A 
throughout the paper. 



II. ANALYTICAL FORMULA FOR 
ESTIMATING COSMOGENIC u INTENSITY 

The neutrino flux per unit energy, dJ v /dE v , is gener- 
ally written as, 
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The first integral represents the total contribution of 
UHECR sources in the redshift up to z max , where z max is 
the maximum redshift of the UHECR source distribution 
or, in other words, the time of the first UHECR emission 
in the Universe. tp( z s) represents the cosmic evolution of 
the spectral emission rate per co-moving volume and hq 
is the number density of UHECR sources at the present 
universe. The relation between time and redshift is given 



| = [H (l+Z a )y/tl M (l + Zs) 3 + OV]- 1 - 



The second integral calculates the total neutrino flux 
expected from a single UHECR source at redshift z s 
generated via UHECR interactions at redshift z u (< z s ). 
dNp^y / dEf,dt 9 is the yield of generated neutrinos with 
energy of E 9 per unit time in the UHECR laboratory 
frame (the CMB rest frame). Suffixes s and g represent 
the quantities at the positions of UHECR sources and 
neutrino generation, respectively. The delta function in- 
dicates the neutrino energy loss due to the expansion of 
the universe. 

The neutrino yield, dN p _> v / ' dE 9 dt 9 , at redshift z v by 
the GZK mechanism is expressed by a convolution of the 
UHECR intensity from a source at z s , the CMB photon 
density, and the photo-pion interaction kinematics as, 
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where dNcn/dEcR is number of the UHECRs per unit 
time and energy at the redshift z v originating from a 
source at z s , and s is the Lorentz- invariant Mandelstam 
variable, the square of invariant mass of the cosmic-ray 
nucleon and the target CMB photon, the photo- 

pion production cross section, dp^^jdE 9 is the energy 
distribution of neutrinos from the photo-produced pion, 
and dn-y/ds is the CMB photon number density in the 
UHECR frame per unit s. 

Wc introduce following approximations to simplify the 
calculations: 1) the contribution of UHECR colliding 
with IR/O background is negligible and only the con- 
tribution of photo-pion production cross-section from A- 
resonance is considered in collisions of UHECRs and 
CMB photons, and 2) the kinematics of the photo-pion 
production is represented by a single pion production. 
The first approximation allows the photon number den- 
sity dn-y/ds to be analytically obtained with the mod- 
ification to the black-body distribution Q. The con- 
tribution of neutrinos induced by UHECR interactions 
with IR/O becomes sizable only in energy region below 
100 PeV [15[ while the effect is small in the higher energy 
region. Similarly the neutrinos from photo-produced pi- 
ons outside the A-resonance are mostly visible only in 
the lower energy range below 100 PeV [Tg, and the sin- 
gle pion production is the most dominant channel in the 
A-resonance. The A-resonance approximation simplifies 
the integral on s in Eq. [3] to a multiplication of the inte- 
grand at s = s_r(~ 1.5 GeV 2 ), where sr is the Lorentz- 
invariant Mandelstam variable at the A-resonance, and 



As#(~ 0.6 GeV 2 ), the width of the A-resonance. The 
2nd approximation then gives [l7| , 
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where — m^/m 2 ~ 0.57 is the muon-to-pion mass 
squared ratio. The factor three arises from the fact that 
three neutrinos are produced from the 7r meson and /i 
lepton decay chain. The allowed range of E 9 due to the 
kinematics is given by, 



E 9 

< -f < 1 
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where neutrino mass is neglected. With a good approx- 
imation that a single pion is isotropically emitted in the 
center-of-momentum frame, one obtains 
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pion normalized by the parent proton energy Ecr- x 
are the maximal and minimal bound of x^ due to the 
kinematics and given by 
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where m p is the proton mass. 

Then we obtain the neutrino yield, Eq. [21 expressed as 
an analytical function with only a single energy integral; 
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Here fee is the Boltzmann constant, T is present temper- 
ature of the CMB. Ea = (sfl —TO 2 )/4fcs? 1 corresponds to 
the energy of UHECR protons colliding via A-resonance 
at the present universe. Suffix R denotes the values at 
the A-resonance in the photo-pion reaction. For example, 
a^ p = 2.1 x 10~ 28 cm 2 represents the photo-pion produc- 
tion cross-section of channel yp —> nir + at the A reso- 
nance and Xft is given by Eq. [6] with s = sr. 

The parameter £r reflects the kinematics bounds, 
Eqs. @]and[ni and is defined by 
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dNcR/dEcn is calculated by the energy loss formula 
with the Continuous Energy Loss (CEL) approxima- 
tion [lU represented by 
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where Agzk is the energy attenuation length governed 
by the GZK mechanism, mainly due to the photo-pion 
production of UHECRs and the CMB. The factor (1 + z) 3 
accounts for the increase of CMB photon number density 
with redshift z. 

Here we introduce the final approximation that the 
energy attenuation length of UHECR by the GZK mech- 
anism, Aqzk, is constant with energies above -Egzk = 
10 20 eV. While Agzk rapidly decreases with cosmic-ray 
energy increase, it turns to be a slight increase or con- 
stant above ~ 3 x 10 20 eV for z v ~ 0. Neutrinos from 
z v > 1 are dominant contribution to the cosmogenic 
neutrino intensity at earth and the turnover energy are 
shifted to lower energies < Eqzk due to the redshift 
effects for the universe of z v > 1 Therefore this 

approximation reasonably describes the UHECR energy 
loss profile to calculate the neutrino yield. Assuming 
the primary UHECR spectrum from a source at z s fol- 
lows the power law described by dNcn(z s ,z s )/dEcR — 
kcr(E C r/ E GZ K)~ a up to E max , the maximal injected 
energy from a source, then the dNcR(z s , z u )/dEcR is an- 
alytically obtained by, 
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where f(z) = + z) 3 + and kcr is a nor- 

malization constant. With Eq. [101 the energy integral 
on Ecr in Eq. [7] becomes an addition of integrals in the 
forms of / dyy-( a+ V h^l-e^) and / dyy~^ a+ ^ ln(l- 
e~ 1 / y )lny. An asymptotic approximation with numeri- 
cal constants is found to provide approximate solutions 
of these integrals. The final formula of the neutrino yield 
is then obtained as, 
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and 

£i = <l xiEAX^(l-r n )< E^(l+z u ) <xiEAX%(l-r w ), 
lO otherwise. 

(13) 

One can also find that xiE&x^l — r w ) in these equa- 
tions represents the effective energy of neutrinos from 
decay of the pions of which energies are kinematically 
allowed energy maximum (x^) and minimum (x^) from 
A-resonance in the 7p collision. E^(l+z u ) = E u (l + z u ) 2 
factor reflects the redshift dependence of the CMB tem- 
perature and the redshift energy loss of neutrinos at z v . 

Eq. [1] with the formula Eq. [TT] finally gives the cosmo- 
genic neutrino flux with double integrals of redshift z s 
and z v . The z v integral is analytically solvable neglect- 
ing 0((AGZK-ffo/ c ) 2 ) or higher order terms as the energy 
attenuation length is much shorter than the cosmological 
time dimension. Then the final form of the cosmogenic 
neutrino intensity is obtained as, 
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Here For represents the UHECR intensity above -Egzk 
and described by 
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assuming that E max ^> Eqzk- Fqr gives the normal- 
ization of the neutrino flux in the present formulation in 
Eq. Q3J It can be estimated by the observational data for 
actual calculation. 

( in Eq. Q3] is the term which account to the redshift 
dependence and given by, 
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where xq — 0.275, and x± = 0.16 are the empirically de- 
termined numerical constants. e arL d £i are either unity 
or null, depending on neutrino energy. These are conse- 
quences of the kinematics bound for pions and neutrinos 
(Eq. [5J , and given by 

eo [l ES(l + Zu)<XiEAX$(l-r v ), (12) 
I otherwise. 



where eo and e% are obtained by Eq. Q2] and [12] respec- 
tively, replacing z v by z s . The cosmic evolution function 
-0(z s ) is now parametrized as (1 + z s ) m such that the 
parameter m represents the "scale" of the cosmic evolu- 
tion often used in the literature [20]. The integral on z s 
in Eq. [TB] is analytically solvable when we use the fact 
that ^/£!m(1 + z s ) 3 + fl\ 3> Q\ in most of the integral 
range. Finally we obtain the final form of the red shift 
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dependent part of the analytical formula £ as, 
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where 7 m = (a + m) — 3/2. z up and Zd OTn are the maxi- 
mum and minimum bounds of the redshifts, respectively. 
These redshift bounds are associated with z max in Eq. [I] 
but also depend on neutrino energies E v due to kinemat- 
ics of 7r-decay and the redshift energy loss. z up is given 
by, 
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Zdown is also given by Eq. [TS] replacing by x R . 

See appendix[B]for the case of the astronomical objects 
of which cosmological evolution become constant above 
a certain redshift (See Ref. [Il[ for example). 



III. VALIDITY OF THE ANALYTICAL 
FORMULA 

The analytical formula for estimating cosmogenic neu- 
trino fluxes (Eqs. Q3J H7J and IT5|) is derived under sev- 
eral assumptions. Here we demonstrate the applicability 
of the formula in estimating neutrino flux in 100 PeV 
< E v < 10 EeV which is the main energy range of sev- 
eral cosmogenic neutrino searches In Table U 
the cosmogenic neutrino integral flux above 1 EeV ob- 
tained by the analytical formula with a — 2.5 are pre- 
sented. We use the UHECR intensity F C r(> E GZ k) — 
2.96 x 10 -21 cm -2 sec -1 sr _1 in the present study, which 
is obtained from the measurement of the HiRes experi- 
ment Q. The fluxes obtained by the full numerical cal- 
culations with the same or comparable source evolution 
parameters are also listed for comparison. The values in 
each parameter subset show an agreement within a fac- 
tor of two for a comparable evolution scenario in the wide 
range of parameter numbers. 

Figure [1] presents the neutrino fluxes obtained with the 
present analytical estimation and the full blown numeri- 
cal calculations. The fluxes calculated with the different 
techniques show the best agreement at ~1 EeV, the cen- 
tral energy in the cosmogenic neutrino search with Ice- 
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Ahlers et al. [9] 
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Ahlers et al. [9] 
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("the best fit") 
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TABLE I. Cosmogenic neutrino fluxes predicted by the 
model-dependent full numerical calculations and those given 
by the present analytical formula with the corresponding pa- 
rameters on source evolution. The numbers by the full calcu- 
lations were converted to be the sum over all three neutrino 
flavors from the original when appropriate. 



Cube [ll| . The present formula provides a reasonable es- 
timate of the neutrino flux from 100 PeV to 10 EeV with 
uncertainty of factor of ^two. Some deviations in the 
analytical formula from the full blown numerical calcu- 
lations arise mainly from the uncertainty in the intensity 
of the extra-galactic UHECR component allowed by the 
observed UHECR spectrum, and the accuracy of the ap- 
proximations used in derivation of the analytical formula. 
We discuss these issues in SccfVl 



IV. RESULTS 

A. The relation between the v flux and the 
cosmological evolution of the sources 

Shown in Fig. [5] is the distribution of the cosmogenic 
neutrino integral fluxes with energies above 1 EeV in 
the parameter space of the evolution of UHECR sources 
(m, z max ) calculated using the derived analytical for- 
mula. The fluxes vary by more than an order of mag- 
nitude with the evolution parameters. The distribution 
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FIG. 1. Integral neutrino fluxes, J [cm -2 sec -1 sr" 1 ], as a 
function of neutrino energy. Bold lines represent the present 
analytical estimates and thin lines represent corresponding 
predictions by the full numerical calculations |3| or the Monte- 
Carlo simulations 0, Q1J 




m 



FIG. 2. Integral neutrino fluxes with energy above 1 EeV, 
J [cm -2 sec -1 sr _1 ], on the plane of the source evolution 
parameters, m and z max - 



demonstrates that the neutrino intensity can indeed be 
an observable to imply the characteristics of the UHECR 
sources. The plot shows that cosmogenic neutrino flux 
around 1 EeV is mostly determined by source emissiv- 
ity history up to redshift of z s ~ 3. This is because the 
contributions of sources at z s > 3 represent only a small 
fraction of the total flux due to the redshift dilution [ll| . 



B. Constraints on UHECR origin with the 
IceCube diffuse neutrino flux limit 

Here we estimate the expected event rates with the 
IceCube neutrino observatory by using the derived an- 
alytical formula. The analytical function is valid in the 
IceCube cosmogenic neutrino detection energy range dis- 
tributed around 1 EeV [11]. Convolution of Eq. [T4l with 
the IceCube neutrino effective area [ill, E3 gi yes the 
event rate for the entire phase space of the evolution pa- 
rameter to and the maximal redshift z max ■ Full mixing in 
the standard neutrino oscillation scenario is assumed and 
the intensity of neutrinos of each of three neutrino flavors 
corresponds to one third of the estimated neutrino inten- 
sity by the analytical function. The Feldman-Cousins 
upper bound [24| then defines the excluded region on the 
m-Zmax plane at a given confidence level. Figure [3] dis- 
plays the resultant constraints. The shaded region rep- 
resents the factor of two uncertainty in the analytical 
estimation discussed in the previous section. The upper- 
limit with the IceCube 2008-2009 data [ll| has already 
started to constraint on hypotheses of UHECR sources 
with strong evolution of m > 4.5. While this bound 
may be still weaker than that by the Fermi diffuse 7-ray 
flux measurement nevertheless the limit by neutri- 
nos is important because the neutrino estimate does not 
involve the uncertainties of the assumptions of E max nor 
the EBL intensity. The full IceCube five-year observation 
would certainly probe the most interesting region of the 
source evolution phase space where the strong candidates 
for the UHECR sources of the powerful astronomical ob- 
jects such as radio galaxies and GRBs are included. 



V. DISCUSSION 

The derived analytical formula to calculate intensities 
of the neutrinos produced by the GZK mechanism in the 
range between 100 PeV and 10 EeV is used to constrain 
the cosmological evolution of the UHECR sources. 

The largest uncertainty in the present analytical for- 
mula at the lower energy range (E v <C 1 EeV) is due to 
the omission of the IR contribution to the cosmogenic 
neutrino production. Photo-produced pions from the 
UHECR interactions with the IR background are major 
origin of neutrinos with energies below 10 PeV, however, 
they are relatively minor in the higher energies where 
we mainly discuss in the aspect of the cosmogenic neu- 
trino detection by IceCube. The amount of the IR con- 
tribution was studied, for example, in the calculations in 
Refs. [IGJl. The study in the former reference exhibits 
much higher contribution of the IR background than in 
the latter where the effect is suppressed mainly due to 
the introduction of the minimal energy of extra-galactic 
UHECR population. These differences can be seen in 
figure [TJ the low energy component in Ref. [HI is sub- 
stantially emphasized compared to the other calculations. 
These variations in the estimation of the IR contribution 
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m m m 

FIG. 3. Constraints on the UHECR source evolution m and z max with the IceCube 2008-2009 flux limit [llj (left) and with 
the full IceCube five-year sensitivity [12] (middle). The area above the solid lines are excluded by null detection of v events. 
The shaded belts represent uncertainties in the present analytical estimation. The right panel shows the full IceCube 5 year 
constraint when the emission rate per co-moving volume becomes constant above z s of 1.0. Excluded region at 68 % confidence 
level (corresponding to ~ 1.1 events assuming the same background rate of the IceCube 2008-2009 measurement [H]]), and 90 
% confidence level (~ 2.2 events) are displayed. 



to the cosmogenic neutrino intensities are considered to 
be an additional uncertainty to the IR background yield 
itself which is also not firmly understood [l9|,|27|. Here we 
would like to emphasize that the omission of the IR back- 
ground leads to a conservative constraint on the UHECR 
source evolution. 

The second largest uncertainty is concerned with -Fcr> 
the UHECR intensity above E GZK ~ 10 20 eV. The works 
in Rcfs. [9|, l28| allowed a sizable variation in the UHECR 
intensity within 99% confidence level of the statistical 
test against the observed data. It indicates that an 
extreme case of the UHECR intensity may lead to a 
large departure from the present estimate of the neutrino 
fluxes. For instance, the difference of their estimate for 
the scenario of (m, z max ) = (2, 2) found in Table fl] arises 
from their assumption of very steep UHECR spectrum 
leading to a minimal -For estimation. This uncertainty 
is, however, expected to be reduced in future when the 
statistical uncertainties in the observations of UHECRs 
and/or the systematic uncertainty in the energy estima- 
tion are improved. 

We would like to also emphasize that the neutrino in- 
tensity below 10 EeV is not largely affected by the de- 
tailed behavior of UHECR proton propagation in extra- 
galactic space. This is because these neutrinos are mostly 
generated at cosmological distances away, which are sub- 
stantially longer than the UHECR proton energy attenu- 
ation length in the CMB field. It is also suggested by no 
explicit dependence of Aqzk in the final formula Eq. 1141 
This is related with the fact that the cosmogenic flux 
below 10 EeV is relatively insensitive to E max and a, 
the maximal injection energy of UHECR protons from 
their sources and the spectral index of UHECR spec- 
trum, respectively, while the flux above 10 EeV is sen- 
sitive to those parameters @, [IH, [23[ . A scan of the pa- 
rameter spaces of the cosmogenic neutrino sources for 
some known classes of astronomical objects with a nu- 



merical Monte- Carlo method was made in Ref. [HI, [l9| 
and it was also shown that the intensity around 1 EeV is 
stable against £ max variation and the transition models 
between the Galactic and extra-galactic cosmic-ray com- 
ponents. These observations are consistent with the fact 
that the neutrino intensities around 1 EeV by the rela- 
tively old works @, [lj| assuming harder UHECR spec- 
trum of a = 2.0 and higher E max , and those by the 
recent works 0, [l5[ with a ~ 2.5 shows an agreement 
also within a factor of two. The difference between the 
present analytical formula and the full blown simulation 
above ~10 EeV in figure [1] is attributed to responses 
to -Emax- The present analytical estimates of neutrino 
fluxes for 100 PeV < E v < 10 EeV, the main energy 
range by the IceCube cosmogenic neutrino search, is ro- 
bust against these parameters. We should note however 
that we use i*cR f°r the normalization constant assum- 
ing i^max Eq-zk- If B m ax is comparable or lower than 
£gzk, the neutrino yield strongly depends on i^max and 
the present simple treatment is not capable of providing 
reasonable estimates of the cosmogenic neutrino fluxes. 

The present analysis indicates that a five-year observa- 
tion by the IceCube observatory will scan the source evo- 
lution parameter space of the most interest where many 
of the proposed UHECR astronomical sources are dis- 
tributed. A null neutrino observation then would imply 
that either UHECR sources are only locally distributed 
{zmax ^ 1), very weakly evolved (m < 3), or the UHE- 
CRs are not proton dominated but heaver nuclei such 
as irons after all. The first two possibilities may lead to 
a speculation about the highest energy particle emission 
from an entirely different and probably dimmer class of 
objects than currently suggested. The last possibility has 
also been discussed with the measurement of the depth of 
maximum of air-showers by the Auger collaboration [2o| . 
Neutrino search in ultra-high energies provides a comple- 
mentary constraint on the proton fraction of UHECRs in 
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this case [2i 



The neutrino flux is then calculated by 



VI. 



SUMMARY 



dt 



(Zs) 



We have derived the analytical formula to estimate the 
cosmogenic neutrino fluxes for wide range of cosmolog- 
ical evolution parameters of UHECR emission sources. 
The analytical formula provides a practical tool for es- 
timating the neutrino intensity at around EeV energy 
region with a limited accuracy within a factor of ~two. 
The obtained analytical estimates have indicated that 
the present IceCube neutrino limit in 100 PeV - 10 EeV 
energies disfavors the scenarios with the strongly-evolved 
UHECR sources. The future IceCube observation will be 
able to scan most of the interesting parameter space of 
UHECR source evolution. Furthermore, while the deep 
and highly energetic part of Universe is inaccessible with 
photons or cosmic-rays due to the CMB field, the cur- 
rent study implies that the neutrinos can be used as a 
rare tool to probe the far Universe. 

With the greater statistics of ultra-high energy neu- 
trino detections by the future neutrino telescopes of 
-100 km 2 areas, such as ARA [H and ARIANNA 
the analytical formula allows us to specify the astro- 
nomical classes of ultra-high energy cosmic-ray sources. 
The pioneer prediction of the cosmogenic neutrinos in 
1960's [4] will finally lead to revealing the characteristics 
of UHECR emission mechanism in a near future. 
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Appendix A: The flux calculation based on 
energetics 

Recently the possible upper-limit on the cosmo- 
genic neutrino flux has been discussed solely using the 
Fermi measurement on extra-galactic diffuse 7-ray back- 
ground [3lj |. In this work, the neutrino flux was approx- 
imately estimated by the energetics argument; calculat- 
ing energy channeling into secondary neutrinos from a 
UHECR proton during its propagation in the CMB field. 



/ dE CR E CR ^R u (E CR )^. (Al) 

J dhj v 

Here dNc R /dEc R is the injected UHECR proton spec- 
trum (~ E~ a ) at the source redshift z s , R v is a fraction 
of UHECR proton injection energy carried by the sec- 
ondary neutrinos, dp v jdE v is a distribution of neutrino 
energy. R v was calculated in the earlier work [l6j repre- 
sented by a numerically fitted function as 



R v = 



0.45 



, ( 2x10" GeV V 

V (1 + 2 S )-E C R ) 



(A2) 



1 



While the original work |3l| represented dp v jdE v as — 
S(E L , — Ecr/ (20(1 + z s ))) approximating each secondary 
neutrino receiving 1/20 of the UHECR proton energy, 
we found that the single pion kinematics approximation 
would give a better agreement in the neutrino spectral 
shape with those obtained by the full blown simulation. 
It is then written as 



dE u 



(l+z s )[E CR (x+^x R )(l-r n )]-Hn(^j , (A3) 



where £r is given by EqfSl replacing Ef, with E v (l + z s ). 

This approach has an advantage that it does not rely 
on the A-resonance approximation. Although we are not 
able to find out a complete analytical solution of the in- 
tegrals in Eq. IA11 the numerical calculations indeed con- 
firmed that the formula Eq . I A 1 I reasonablv reproduces the 
full simulation/numerical calculation results. It gives a 
better agreement than our formula at around 100 PeV, 
owing to inclusion of direct pion production yielding a 
pair of 7r + 7r~ by Eq. IA2I [lfj. However, this energetics- 
based formulation significantly overestimates neutrino in- 
tensities with energy above 1 EeV. We suspect it due to 
the neglecting of the energy loss of UHECR protons. En- 
ergy of an UHECR proton is in many cases largely lower 
than its injected energy when it yields neutrinos, because 
of energy loss by the photo-pion production during the 
UHECR propagation. Without accounting this effect, 
higher energy neutrino production is over weighted in 
the formulation. The overproduced high energy neutri- 
nos are then red-shifted and accumulated even in PeV 
regime when UHECR sources are strongly evolved. As 
a consequence, the estimated intensity departs from the 
calculation with the full blown simulation in case of the 
strong evolution scenario. Since the GZK neutrino search 
by the IceCube is sensitive to EeV range and emission 
from strongly evolved sources, we concluded that the 
energetics-based formulation is not appropriate for our 
purpose. 



Appendix B: The analytical formula for the partially 
constant source evolution 



Some classes of astronomical objects, like the galaxy 
star formation, seem to exhibit evolution nearly constant 
above a certain redshift. For these cases, cosmological 
evolution of sources are written as, ip(z s ) ~ (1 + z s ) m for 
< z s < z and ip(z s ) ~ (1 + z) m for z < z s < z max . The 
redshift dependence term ( (Eq. flBI) is then obtained with 
minor modifications on the formula Eq. 1171 and given as, 



tions as Eq. [TTlas. 

C = e" 2 — n-*(i + *) m 

(!Jm(1 + ~Zu P f + Ov)^ (x- (a+1) In ($) + 



-(«+3) -^r 



In I xi 



E u (l + z up ) 



+- 



7« 



— (Am(1 + Zdown) + Oa) 3 X 



(c+1) 



hi 



-(njitf(l + 2 doMn ) 3 + nA)^asr (a+3) e »i 



In I a;i 



Ea 



E v (l + Zdou 



r^ji(l -*v) ) + 



7a 



(B2) 



C = C* + C 



(Bl) 



where £j is given by Eq. [17] with replacing z mQX by z, 
accounting for the evolution up to z s < z. 



The additional term £ is obtained in the similar func- 



by 



Where j a = a — § and the redshift bounds z up is given 



1 + z u 



1 + 5 



55 S 



' i% 4( 1 -^)<^< 



(1+4)2 



1 H - 2>max 



E v < 



"TT+lp" 

1 i E A 1 a('- r ») 



(B3) 



Zdown is written as the same Eq. IB3I with replacing 
by Xr. 
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